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Helicity transfer in a shell model of turbulence is inves- 
tigated. We show that a Reynolds-independent helicity flux 
is present in the model when the large scale forcing breaks 
inversion symmetry. 

The equivalent in Shell Models of the "2/15 law", obtained 
from helicity conservation in Navier-Stokes eqs., is derived 
and tested. 

The odd part of helicity flux statistic is found to be domi- 
nated by a few very intense events. 

In a particular model, we calculate analytically leading and 
sub-leading contribution to the scaling of triple velocity cor- 
relation. 

One of the most intriguing problems in three dimen- 
sional fully developed turbulence (FDT) is related to the 
appearance of anomalous scaling laws at high Reynolds 
numbers, i.e. in the limit when Navier-Stokes dynamics 
is dominated by the non-linear interactions. 

The celebrated 1941 Kolmogorov theory (K41) was 
able to capture the main phcnomenological ideas by per- 
forming dimensional analysis based on the energy trans- 
fer mechanism. Kolmogorov postulated that the energy 
cascade should follow a self-similar and homogeneous 
process entirely dependent on the energy transfer rate, 
e. This idea, plus the assumption of local isotropy and 
universality at small scales, led to a precise prediction 
on the statistical properties of the increments of turbu- 
lent velocity fields: 5v(l) ~ \v(x + Z) - v{x)\ ~ (Ze(Z)) 1 / 3 . 
From this the scaling of moments of Sv(l), the structure 
functions, can be determined in terms of the statistics of 
e(Z), i.e. 



S p (l) = ((M0) P > = C P Ml))^)P/ 3 



(1) 



where C p are constants and the scale I is supposed to be 
in the inertial range, i.e. much smaller than the integral 
scale and much larger than the viscous dissipation cutoff. 
If S P {1) ~ JCCp) and (<*(/)) ~ Z T ^ then 

C(p)=p/3 + r(p/3). 



(2) 



In the K41 the e(Z) statistic is assumed to be Z- 
independent, or r(p) = 0, implying = |, Vp. On the 
other hand, there are many experimental and numerical 
|fl,pl results telling us that K41 scenario for homogeneous 
and isotropic turbulence is quantitatively wrong. Strong 



intermittent bursts in the energy transfer have been ob- 
served and non trivial r(p) set of exponents measured. 

Moreover, the problem of investigating scaling proper- 
ties of observables with the same physical dimensions but 
with different tensorial structures has been only recently 
addressed |J|. 

Many different authors have focused their attention 
on the possible role played by helicity, the second global 
invariant of 3d Navier-Stokes eqs. [^|-|8), for determining 
leading or sub-leading scaling properties of correlation 
functions in the inertial range. 

Recently [0J1, an exact scaling equation for the third 
order velocity correlations entering in the helicity flux 
definition has been derived under two hypothesis: (i) 
there exists a non-vanishing helicity flux, (ii) the flux be- 
comes Reynolds independent in the limit of FDT. This 
relation predicts a r 2 scaling for a particular third order 
velocity correlation. The new relation has been called 
"2/15 law" because of the coefficient appearing in front 
of the r 2 in analogy with the "4/5 law" derived by Kol- 
mogorov for the third-order structure functions entering 
in the expression of energy flux. In the "4/5" law the 
scaling of a different third-order velocity correlation is 
found to be linear in r. 

This simple fact tells us that different velocity correla- 
tion with the same physical dimension but with different 
tensorial structure may have very different scaling prop- 
erties. 

Moreover, even if overwhelming evidences indicate that 
the main physics is driven by the energy transfer, there 
can be some sub-leading new intermittent statistics hid- 
den in the helicity flux properties. 

Homogeneous and isotropic turbulence has, by defini- 
tion, always a vanishing mean helical flux. Nevertheless, 
both fluctuations about the zero-mean in isotropic cases 
and/or net non-zero fluxes in cases where inversion sym- 
metry is explicitly broken can be of some interest for the 
understanding of fully developed turbulence. 

In this letter, we analyse the helical transfer mecha- 
nism in dynamical models of turbulence |lO|Jq,pT[|, built 
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such as to explicitly consider helicity conservation in the 
inviscid limit. 

We give the first strong numerical evidence that a 
Reynolds-independent helicity flux is present in cases 
where the forcing mechanism explicitly breaks inversion 
symmetry. 

We confirm that in all cases where two fluxes can 
coexist in the inertial range, velocity correlations with 
the same physical dimension but with different transfor- 
mation properties under inversion symmetry, can show 
strongly different scaling behaviour. 

In the following, we briefly summarize the main moti- 
vation behind the introduction of Shell Models for turbu- 
lence. We present the equivalent derivation of the "2/15 
law" for the helicity flux in our Shell Model language 
and we conclude with our numerical results about the 
Reynolds independency of helicity flux and about its sta- 
tistical intermittent properties. 

Shell models have demonstrated to be very useful for 
the understanding of many properties connected to tur- 
bulent flows jjJI- The most popular shell model, the 
Gledzer-Ohkitani-Yamada (GOY) model ( @), has 
been shown to predict scaling properties for (for a 
suitable choice of the free parameters) similar to what is 
found experimentally. 

The GOY model can be seen as a severe truncation of 
the Navier-Stokes equations: it retains only one complex 
mode representative of all Fourier modes in the 

shell of wave numbers k between k n = ko2 n and k n +\. 

It has been pointed out that GOY model conserves 
in the inviscid, unforced limit two quadratic quantities. 
The first quantity is the energy, while the second is the 
equivalent of helicity in 3D turbulence Jl7[ |. In two re- 
cent works |8|,|ll| the GOY model has been generalized 
in terms of shell variable, it+,w~, transporting positive 
and negative helicity, respectively. It is easy to realize 
that only 4 independent classes of models can be derived 
such as to preserve the same helical structure of Navier- 
Stokes equations @. All these models have at least one 
inviscid invariants non-positive defined which is similar 
to the 3D Navier-Stokes helicity. In the following, we will 
focus on the intermittent properties of a mixture of two 
of such a models. The mixture is a linear combination of 
the old GOY model (extended to have u + ,u~) plus an- 
other model which has a different helical interaction and 
which has already been extensively investigated |0,E9[. 
We focus only on two of the four possible models because 
they are the only two classes of models which show a clear 
forward energy cascade (see for more details). The 
time evolution for positive-helicity shells reads JO]: 



u+ = ik n (A n [u, u] + xB n [u, u])* — vk n u\ + 



3n,n f j 



(3) 



with the equivalent eqs, but with all helicity signs re- 
versed, for ii~ . In (H), x defines the relative weights of 



the two models in the mixture, v is the molecular vis- 
cosity, /+, / _ are the large scale forcing and A[u, u] and 
B[u,u] refer to the non- linear terms of the two models. 
Namely: 

A n [u, u] = u~ +2 u+ +1 + b 3 u^ +1 u+_ 1 + c 3 u-_ lU -_ 2 , (4) 
B n [u, u] = u+ +2 u n+1 + &iu^ +1 u^_ a + c 1 u; i _ 1 u+_ 2 . (5) 

It is easy to verify that for the choice b 3 = —5/12, c 3 — 
— 1/24, &i = — 1/4, ci = —1/8 there are only two global 
inviscid invariants p| : the energy, E 
|u"| 2 )), and helicity, H = J2 n Ml"™ ? ~ \ u r, 

The equations for the fluxes throughout shell number 
n are: 



I 2 ) 
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d/dt J2 E i = k n < ( 



llliu)„ > 



Ei. 



(6) 



i—l,n 



d/dt Hi = kl < (uuu)* > -v Y, H > + H ™ ( 7 ) 



where Ei and Hi are the energy and helicity of the nth 
shell, respectively: E t —< \u+\ 2 + \u^\ 2 >, H, t = ki(< 
\ u t\ 2 — l u rl 2 >)■ and H in are the input of energy 
and helicity due to forcing effects, E in = K(< / + (u^)* + 

/-(«r)* >)> H in = < /+«)* - /"(%)* >)■ 

In (||) and (0) we have introduced the triple correlation: 

< > = (A+ +1 + A~ +1 ) + (h + 1/2)(A+ + A") 

+ x ((r+ +1 + r- +1 ) + (&! + i/2)(r+ + r-)) , (8) 

< (uuu)* > - (A+ +1 - A- +1 ) + (6 3 + 1/4) (A+ - A-) 

+ x ((r+ +1 - r- +1 ) + (h + 1/4) (r+ -r~)) , (9) 

and 



A,| = ($s(u+ +1 u n , T+ = (S(u n+2 u+ +1 u+)) . 

(10) 

Assuming that there exists a stationary state we have 
= = 0, where Ilf = k n ((uuu)%) and 

= k n ((uuu)^). Moreover, in the inertial range we 
can neglect the viscous contribution in (@) and ([7]), ob- 
taining: 

< (uuu)% >= k~ l E in , < (uuu)n >= K 2 H m . (11) 

Supposing that there exist the energy and helicity fluxes, 
Ei n = e, Hi n — h (the latter different from zero only 
if / + 7^ / — ) and supposing that both are Reynolds- 
independent, we have in the inertial range: 



< (uuu)^ >^ fc n 1 , < {uuu)n >^ k n 



(12) 



Relation ( jll| ) is the equivalent of what found for helical 
Navier-Stokes turbulence in |||]]). 

Figure 1 reports results for the helicity flux in a numer- 
ical simulations done with two different Reynolds num- 
bers, Re ~ 10 5 , Re ~ 10 9 for a choice of mixture parame- 
ter x = 0.1 and other numerical inputs as follows, N = 16 
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and 26, v = 10 5 and 2-10 9 . A non-zero helical flux was 
obtained using a forcing term breaking inversion symme- 
try: /+ = 5 • 1CT 3 (1 + i), f - = /+/10. Time marching 
was obtained by using a slaved Adams-Bashforth algo- 
rithm for a number of iterations equal to several thou- 
sands of the typical eddy turnover time. A clear inertial 
range with a non- vanishing helicity flux is detected. The 
extension of the range where helicity flux is roughly con- 
stant scales with the Reynolds number. Moreover, the 
flux intensity is roughly constant at changing Reynolds 
number, giving the first evidence that the model can si- 
multaneously support both energy and helicity transfers 
and that both of them are Reynolds- independent. 

Let us remark that this is only possible due to the 
non-positiveness of helicity; in 2d turbulence, for exam- 
ple, a similar results, concerning enstrophy and energy 
cascades, is clearly apriori forbidden. 

As for the statistics of helicity transfer, we measured 
the scaling exponents of the moments of energy and he- 
licity fluxes: 

4 P) ^< (W?) W3) >~ (is) 

E|)=<((nnO^/ 3 )>~fc-^). (14) 

As one can see in figure 2 we have that the even part of 
the two statistics coincides, i.e. ((2p) = ip(2p). On the 
other hand, the scaling exponents of odd moments are 
different. 

The difference in odd moments is the signature of 
strong cancellation effects in the statistics connecting 
fluctuations at different scales. The picture we have in 
mind is that the main effect driving turbulent fluctua- 
tions is due to the energy cascade process, with its in- 
termittent fluctuations measured by the £(p) exponents. 
Superimposed to the energy transfer, there are "topolog- 
ical" fluctuations introduced by the asymmetric forcing 
and measured by the odd-part of the helicity turbulent 
transfer. 

Let us notice that helicity flux fluctuations are much 
larger then the average helicity flux. This clearly distin- 
guishes the helicity transfer mechanism from the energy 
transfer mechanism. 

In the former, the strong intermittent behaviour shown 
by odd moments tell us that the odd part of the statistic 
is dominated by a few very singular structures. 

In the case of no-mixture (x — 0) one can also ex- 
actly calculate sub-leading scaling for the triple correla- 
tion 5+ = k n (A+ +A~), 6- = kl(A+ -A"). Indeed, 
from expressions (|^) one obtains after some simple alge- 
bra: 



St 



2e- 



1 



V 



n+l 



Ah- 



1 



7 n+l 



1 - Z 



(15) 



where y = -2(b 3 + 1/2) and z = -4(6 3 + 1/4). Being 
both y and z with modulus less then one we recover the 



asymptotic scaling (|12j) and one can also control sub- 
leading correction to it: 
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y" 



1 - 2 ,-i 



y n+l 



-) + k-\2h- 



and 



n+l 



K \2h 



1 - z' 



z 



1-y ' » * 1 

In conclusion, we have studied a Helical shell model 
for turbulence with a forcing which explicitly breaks in- 
version symmetry at large scale. 

We have checked that a Reynolds-independent helicity 
flux establishes in the system, giving the first evidence 
of very different scaling for triple correlations entering in 
the energy flux and helicity flux definitions. 
The odd part of the helicity flux statistic is found to be 
strongly intermittent. 

For a particular class of models we can also calculate ex- 
plicitly sub-leading corrections to pure scaling behaviour 
of typical triple correlation functions. 
The existence of sub-leading terms explicitly tell us that 
scaling laws in turbulent flows must be studied on cor- 
relation functions which have a pure projection on the 
physical relevant quantities. 

There are not reasons why similar effects should not be 
present also in true Navier-Stokes eqs. For example, spu- 
rious intermittent corrections can be detected in cases 
where isotropy is globally or locally violated (as in bound- 
ary layers). 
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FIGURE CAPTIONS 

• Figure 1: 

Helicity-flux, (k^ < (uuu)^ >), versus k n for N — 
16, v = 1(T 5 (dashed line) and N = 26, v = 2-10" 9 
(continuos line). 

• Figure 2: 

Anomalous exponents for the helicity flux, tjj p (cir- 
cles), and for the energy flux, ( p (squares), obtained 
for N = 26 and v = 2 • 1CT 9 . 
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